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Abstract: The Standard Model gauge group can be extended with minimal matter con¬ 
tent by introducing anomaly free U(l) symmetry, such as Lg — L^ or Le — Lr- If the neutral 
gauge boson corresponding to this abelian symmetry is ultra-light, then it will give rise to 
flavor-dependent long-range leptonic force, which can have significant impact on neutrino 
oscillations. For an instance, the electrons inside the Sun can generate a flavor-dependent 
long-range potential at the Earth surface, which can suppress the —)■ Ug appearance prob¬ 
ability in terrestrial experiments. The sign of this potential is opposite for anti-neutrinos, 
and affects the oscillations of (anti-)neutrinos in different fashion. This feature invokes fake 
CP-asymmetry like the SM matter effect and can severely affect the leptonic CP-violation 
searches in long-baseline experiments. In this paper, we study in detail the possible im¬ 
pacts of these long-range flavor-diagonal neutral current interactions due to Le — sym¬ 
metry, when (anti-)neutrinos travel from Fermilab to Homestake (1300 km) and CERN to 
Pyhasalmi (2290 km) in the context of future high-precision superbeam facilities, DUNE 
and LBNO respectively. If there is no signal of long-range force, DUNE (LBNO) can place 
stringent constraint on the effective gauge coupling ae^ < 1.9 x 10“^^ (7.8 x 10“®^) at 
90% C.L., which is almost 30 (70) times better than the existing bound from the Super- 
Kamiokande experiment. We also observe that if > 2 x 10“^^, the CP-violation discov¬ 
ery reach of these future facilities vanishes completely. The mass hierarchy measurement 
remains robust in DUNE (LBNO) if aefi < 5 x 10“^^ (10“^^). 
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1 Introduction and Motivation 

Over the past two decades or so, active attempts have been made both in experimental 
and theoretical fronts to improve our knowledge about neutrinos [1-3]. The three most 
important fundamental issues that have taken center stage in neutrino physics as a part 
of these activities are the following. First issue: how tiny is the neutrino mass? Second 
issue: can a neutrino turn into its own anti-particle? Third issue: do different neutrino 
flavors oscillate into one another? To shed light on the first issue, recently, the Planck 
Collaboration has reported an upper limit on the sum of all the neutrino mass eigenvalues 
of < 0.23 eV at 95% C.L. in combination with the WMAP polarization and baryon 

acoustic oscillation (BAO) measurements [4]. Here, the sum runs over all the neutrino mass 
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eigenstates which are in thermal equilibrium in the early Universe. As far as the second 
issue is concerned, the hunt for neutrinoless double beta decay process is still on which 
violates the total lepton number and requires Majorana neutrinos [5-8]. The third question 
has been answered positively only in recent years [9-11] and now, we have compelling 
evidence in favor of neutrino flavor oscillation [12, 13], suggesting that leptonic flavors are 
not symmetries of Nature. This entails that neutrinos are massive and they mix with each 
other, providing an evidence for physics beyond the Standard Model (SM) [14, 15]. This 
milestone has been achieved due to several fantastic world-class oscillation experiments 
involving neutrinos from the Sun [16-22], the Earth’s atmosphere [23-25], nuclear reactors 
[26-36], and accelerators [37-46] which have enabled us to obtain a clear understanding of 
the lepton mixing pattern in three-flavor scenario. Using the standard Particle Data Group 
convention [14], we parametrize the vacuum Pontecorvo-Maki-Nakagawa-Sakata (PMNS) 
matrix [12, 47, 48] in terms of the three mixing angles: 0 i 2 , O 23 , ^ 13 , and one Dirac-type 
CP phase (5cp (ignoring Majorana phases). In a three-flavor framework, the transition 
probability also depends on two independent mass-squared differences: — rnf 

in the solar sector and Am ‘^2 = ~ ^2 ™ atmospheric sector where m 3 corresponds 

to the neutrino mass eigenstate with the smallest electron component. The smallest lepton 
mixing angle ^13 connects these two sectors and determines the impact of the sub-leading 
three-flavor effects [49-51]. All the neutrino oscillation data available till date have been 
explained quite successfully in terms of these mass-mixing parameters [52-54], excluding 
few anomalous results obtained at very-short-baseline experiments [55]. 

Neutrinos acquire additional phases while travelling in matter, the so-called ‘MSW 
effect’ [56-59] which determined the ordering of the 1-2 mass splitting using solar neutrinos. 
In light of large @ 13 , matter effects are also going to play an important role in presently 
running and future long-baseline [49, 51, 60] neutrino oscillation experiments to settle the 
remaining unsolved issues, namely, the neutrino mass hierarchy^, possibility of leptonic 
CP-violation if (5cp differs from both 0° and 180°, and the octant degeneracy of 623 [61] 
if 023 turns out to be non-maximal. The combined data from the current off-axis 
appearance experiments, T2K [62, 63] and NOpA [64-66], hold promise of providing a 
first hint for these missing links for only favorable ranges of oscillation parameters [67-73]. 
Hence, future facilities consisting of intense, high power, on-axis wide-band beams and large 
smart detectors have been proposed to cover the entire parameter space at unprecedented 
conhdence level [60]. Future superbeam long-baseline facilities with liquid argon detectors, 
Deep Underground Neutrino Experiment (DUNE) [74-78] in the United States with a 
baseline of 1300 km from Fermilab to Homestake mine in South Dakota and Long-Baseline 
Neutrino Oscillation Experiment (LBNO) in Europe involving a path length of 2290 km 
between CERN and Pyhasalmi mine in Finland [79-82] are the two major candidates in 
this roadmap which are capable enough to claim the discovery for the above mentioned 
issues [50, 51]. For both the DUNE and LBNO baselines, the matter effects are quite 
significant which break the eight-fold degeneracies [83, 84] among the various oscillation 


^Two patterns of neutrino masses are possible: m3 > m2 > mi, called normal hierarchy (NH) where 
Am|i > 0 and m2 > mi > m3, called inverted hierarchy (IH) where Am|i < 0 . 
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parameters and improve the physics reach by considerable amount. 

Apart from the SM Vh-exchange interaction in matter, there may well be flavor- 
dependent, vector-like, leptonic long-range force (LRF), like those mediated by the Lg— 
gauge boson which is very light and neutral, leading to new non-universal flavor-diagonal 
neutral current (FDNC) interactions of the neutrinos which can give rise to non-trivial 
three-neutrino mixing effects in terrestrial experiments, and could affect the neutrino prop¬ 
agation through matter [85-89]. Can we constrain/discover these long-range FDNC inter¬ 
actions in upcoming long-baseline neutrino experiments? If this LRF exists in Nature, can 
it become fatal in our attempts to resolve the remaining unknowns in neutrino oscilla¬ 
tion? In this paper, we attempt to address these pressing issues in the context of future 
high-precision superbeam facilities, DUNE and LBNO. 

This paper is organized as follows. We start section 2 with a discussion on flavor- 
dependent LRF and how it arises from abelian gauged Lg — symmetry. Then, we 
discuss the strength of the long-range potential Ug^/gr at the Earth surface generated by 
the electrons inside the Sun. This is followed by a brief discussion on the current constraints 
that we have on the effective gauge couplings Og^^ er of the Lg — symmetry from various 
neutrino oscillation experiments. In section 3, we study in detail the three-flavor oscillation 
picture in presence of long-range potential. We derive the compact analytical expressions 
for the effective oscillation parameters in case of the Lg — symmetry, and also present a 
simple expression for the resonance energy, where 1-3 mixing angle in matter becomes 45° 
in the presence of long-range potential. Next, we demonstrate the accuracy of our approx¬ 
imate analytical probability expressions by comparing it with the exact numerical results. 
We also discuss some salient features of the appearance and disappearance probabilities 
(calculated numerically) for the Fermilab-Homestake and CERN-Pyhasalmi baselines in 
the presence of long-range potential towards the end of this section. We give the similar 
plots for the anti-neutrino case in appendix A. At the beginning of section 4, we give a brief 
description of the main experimental features of the DUNE and LBNO set-ups. Then, we 
study the impact of the long-range potential due to Lg — L^ symmetry on the expected 
event spectra and total event rates for the DUNE and LBNO experiments. In section 5, 
we describe our simulation strategy. Next, we derive the expected constraints on ag^ from 
DUNE and LBNO in section 6.1, and estimate the discovery reach for in section 6.2. In 
section 6.3, we study how the long-range potential affects the CP-violation search of these 
future facilities. Section 6.4 is devoted to assess the impact of LRF on mass hierarchy 
measurements. Finally, we summarize and draw our conclusions in section 7. 

2 Flavor-Dependent Long-Range Forces from Gauged U{1) Symmetries 

The SM gauge group SU{3)c x SU{2 )l x U{1)y can be extended with minimal matter 
content by introducing anomaly free U{1) symmetries under which the SM remains invari¬ 
ant and renormalizable [90]. There are three lepton flavor combinations: Lg — L^, Lg — L,-, 
and L^ — Lj- which can be gauged in an anomaly free way with the particle content of the 
SM [91-93]. These U{1) gauge symmetries have to be broken in Nature in order to allow 
the different neutrino species to mix among each other giving rise to neutrino oscillation 
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as demanded by the recent data [14, 52-54]. This can be understood from the following 
example. If we assume that neutrino masses are generated by effective five-dimensional op¬ 
erator following say, seesaw mechanism, then this operator would remain invariant under 
these U{1) symmetries if they are exact. This will ultimately give us an effective neutrino 
mass matrix containing a Dirac and a Majorana neutrino which remain unmixed and there 
will be no neutrino oscillation. For more discussions on these issues, see references [85, 88]. 
Now, there are two possibilities for the extra gauge boson associated with this abelian 
symmetry^: either it can be very heavy or it is very light but in both the cases, it couples 
to matter very feebly to escape direct detection. If the electrically neutral gauge boson 
{Z') corresponding to a flavor-diagonal generator of this new gauge group is relatively 
heavy^, then the phenomenological consequences can be quite interesting [90-93, 100]. On 
the other hand, if the mass of the gauge boson is very light, then it can introduce LRF 
with terrestrial range (greater than or equal to the Sun-Earth distance) without call upon 
extremely low mass scales [85, 86]. This LRF is composition dependent (depends on the 
leptonic content and the mass of an object) and violate the universality of free fall which 
could be tested in the classic lunar ranging [lOI, 102] and Eotvos type gravity experiments 
[103, 104]. In fact, this idea was given long back by Lee and Yang [105] and later, using 
this idea, Okun [106, 107] gave a 2a bound of a < 3.4 x 10“^® {a denotes the strength of 
the long-range potential) for a range of the Sun-Earth distance or more. Now, the light 
gauge boson Z' should have a mass mz' ~ g {v) where g is the gauge coupling of the new 
U{1) symmetry and it should be < 0.6 x 10“^^ since g ~ ^/a and {v) is the vacuum ex¬ 
pectation value of the symmetry breaking scale. If the range of the LRF is comparable to 
the Sun-Earth distance (~ 1.5 x 10^^ cm) then mz' ~ 1/(1.5 X 10^3 ^ 13 X 10-18 eV 

which means {v) > 2 MeV. 

2.1 Abelian Gauged Lg — ^ Symmetries 

If the extra 17(1) corresponds to Lg — or Lg — flavor combination^, the coupling of the 
solar electron to the Lg—Lfj,^r gauge boson generates a flavor-dependent long-range potential 
for neutrinos [108-110] and can have significant impact on neutrino oscillations [85-89] 
inspite of such stringent constraints on a. This is caused due to the fact that the (Lg—L^, t)- 
charge of the electron neutrino is opposite to that of muon or tau flavor, leading to new 
non-universal FDNC interactions of the neutrinos in matter on top of the SM VF-exchange 
interactions between the matter electrons and the propagating electron neutrinos. These 
new flavor-dependent FDNC interactions alter the ‘running’ of the oscillation parameters 
in matter by considerable amount [111]. Another important point is that the large number 
of electrons inside the Sun and the long-range nature of interaction balance the smallness 
of coupling and generate noticeable potential. As an example, the electrons inside the Sun 

^Models with these symmetries necessarily possess a Higgs sector which also discriminates among dif¬ 
ferent lepton families [94], but we will only focus here on the effect of the extra gauge boson. 

®For an example, the Z' in gauged Le — can be produced in e''"e“ collisions and subsequently decay 
into e^e” or or pairs and can be constrained using the data from LEP/LEP2 [95-99]. 

^Due to the absence of muons or tau leptons inside the Sun or Earth, we do not consider gauged — Lt 
symmetry here. 
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give rise to a potential V'e^/gr at the Earth surface which is given by [85, 88] 

Ve ^/ eriRsE ) = ~ 1-3 X IQ-l'eV , (2.1) 

/ 

where aegleT — deg/er g^uge Coupling of the Lg — symmetry. In 

Eq. (2.1), iVe (~ 10®^) is the total number of electrons inside the Sun [112] and Rse is the 
Sun-Earth distance ss 1.5 x 10^^ cm = 7.6 x 10^® GeV“^. Here af.^jf.T can be identified 
as the ‘fine-structure constant’ of the tl(l)Le-i^.,r symmetry and its value is positive^. 
The corresponding potential due to the electrons inside the Earth with a long-range force 
having the Earth-radius range {Re ^ 6400 km) is roughly one order of magnitude smaller 
compared to the solar long-range potential and can be safely neglected® [85, 88]. The 
long-range potential Eg/i/er iii Eq. (2.1) appears with a negative sign for anti-neutrinos and 
can affect the neutrino and anti-neutrino oscillation probabilities in different fashion. This 
feature can invoke fake CP-asymmetry like the SM matter effect and can influence the 
CP-violation search in long-baseline experiments. This is one of the important findings of 
our paper and we will discuss this issue in detail in the later section. Now, it would be quite 
interesting to compare the strength of the potential given in Eq. (2.1) with the quantity 
Am?/2E which governs the neutrino oscillation probability. Eor long-baseline neutrinos, 
Am^/2E rsj 10 eV (assuming Am^ ~ 2 X 10 ^ eV^ and E rsj 1 GeV) which is comparable 
to Vef,/er even for Uef^/er rsj 10 and can affect the long-baseline experiments significantly 
which we are going to explore in this paper in the context of upcoming facilities DUNE 
and LBNO. 

2.2 Existing Phenomenological Constraints on Lg — Parameters 

There are phenomenological bounds on the effective gauge coupling of the Lg — 

abelian symmetry^ using data from various neutrino oscillation experiments. It was shown 
in [85] that Lg — potential at the Earth due to the huge number of electrons inside the 
Sun suppresses the atmospheric neutrino —)• Vr oscillations which enabled them to place 

tight constraints on using the oscillation of multi-GeV neutrinos observed at the 

Super-Kamiokande (SK) experiment. They obtained an upper bound of Og^ < 5.5 x 10“®^ 
and agr < 6.4 x 10“®^ at 90% C.L. [85]. In [88], the authors performed a global fit to the 
solar neutrino and KamLAND data including the flavor-dependent LRE. They quoted an 
upper bound of ag^ < 3.4 x 10“®^ and < 2.5 x 10“®^ at Su G.L. assuming 6*13 = 0° 
[88]. A similar analysis was performed in [87] to place the constraints on LRE mediated by 
vector and non-vector (scalar or tensor) neutral bosons where the authors assumed one mass 
scale dominance. The proposed 50 kt magnetized iron calorimeter (ICAL) detector at the 

®In our work, we consider the case of a light vector boson exchange which makes sure that ae/j/er is 
positive. It means that for an example, the force between an isolated electron and Vf. is repulsive. 

®The possibility of the local screening of the leptonic force (generated due to the solar electrons) by the 
cosmic anti-neutrinos is also negligible over the Sun-Earth distance [85]. 

^Flavor-dependent long-range leptonic forces can also be generated via the unavoidable mixing of light 
Z' boson of the — Lt symmetry with the Z boson of the SM. This issue was discussed in the context of 
the MINOS long-baseline experiment in [113, 114). 
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India-based Neutrino Observatory (INO) can also probe the existence of LRF by observing 
the atmospheric neutrinos and anti-neutrinos separately over a wide range of energies 
and baselines [89]. With an exposure of one Mton.yr and using the muon momentum as 
observable, ICAL would be able to constrain ^ 1-65 x 10“^^ at 3 ct C.L. [89]. 

3 Three-Flavor Oscillation Picture in Presence of Long-Range Potential 

In this section, we focus our attention to study the impact of flavor-dependent long-range 
leptonic potential (due to the electrons inside the Sun) at the Earth surface when neutrinos 
travel through the Earth matter. In a three-flavor framework, the long-range potential 
of Eq. (2.1) due to Lg — symmetry modifies the effective Hamiltonian for neutrino 
propagation in Earth matter in the flavor basis to 
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where U is the vacuum PMNS matrix [12, 47, 48] which can be parametrized as 

u = i?23(^*23,0) Ri 3(^13, (^cp) -^12(012,0) . (3.2) 

In the above equation, E is the neutrino energy and Vcc is the Earth matter potential 
which appears in the form 

Fee = V2GFNe ~ 7.6 X Eg X „ eV, (3.3) 

lO^^g/cm'* 

where Gp is the Fermi coupling constant, Ng is the electron number density inside the 
Earth, p is the matter density, and Eg = is the relative electron number density. Np, 

Nn are the proton and neutron densities in Earth matter respectively. In an electrically 
neutral, isoscalar medium, we have Ng = Np = Nn and Yf, comes out to be 0.5. In 
Eq. (3.1), Vefi is the long-range potential due to Le — Lp symmetry and its strength is given 
by Eq. (2.1). In case of Le — Lr symmetry, the contribution due to long-range potential 
in Eq. (3.1) takes the form Diag (Igr, 0, —14,-). Note that the strength of the long-range 
potential Vepier does not depend on the Earth matter density® and takes the same value 
for any baseline on the Earth. In case of anti-neutrino propagation, we have to reverse the 
sign of Vcc-, Vep (or Ver), and the CP phase (5cp. 

To judge the impact of the long-range potential Vep in long-baseline experiments with 
multi-GeV neutrinos, we need to compare its strength with the two other main terms in 
Eq. (3.1) which are 2 e^ the Earth matter potential Vcc- hr Table 1, we compare the 
strengths of these three quantities which control the ‘running’ of the oscillation parameters 
in matter. The first oscillation maximum for the DUNE (LBNO) experiment occurs at 
2.56 GeV (4.54 GeV) assuming = 2.44 x 10“® eV^ (see Table 2). In the fourth 

®This is in contrast to the usual non-standard interactions (NSI’s) whose strengths are proportional to 
the Earth matter density. 
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Set-up 

1®* osc. max. (GeV) 

^(eV) 

Ucc (eV) 

Ve^, (eV) 

ae^ = 10 32 

= 10-33 

DUNE 

2.56 

4.8x10-13 

1.1x10-13 

1.3x10-13 

1.3x10-1^1 

LBNO 

4.54 

2.7x10-13 

1.3x10-13 

1.3x10-13 

1.3x10-14 


Table 1. The second column shows the first oscillation maxima for Fermilab-Homestake and 
CERN-Pyhasalmi baselines considering = 2.44 x 10“^ eV^. The fourth column depicts 

the strength of the Earth matter potentials for these baselines. The estimate of the long-range 
potentials 14/i for two different values of is given in the last column. 


Parameter 

Best-fit and In error 

True value 

Marginalization range 

sin^ 012 

0 304'’'® ®13 

U.OUU_g Q32 

0.304 

Not marginalized 

sin^ 013 

n 091 c-to.ooio 

0.02i8_Q gQ3Q 

0.0218 

Not marginalized 

sin^ 023 

U.UOZ_g Q 28 

0.50 

[0.38, 0.64] 

<^cp/° 

306lf^ 

[0, 360] 

[0, 360] 

A?n2i 

10-5eV^ 

'j (r+0.19 
' *^-0.17 

7.50 

Not marginalized 

Amgi 

10-3eV^ 

z.uor_g Q47 

2.44 

Not marginalized 

10-3eV^ 

2 410'''®'®31 

Z.UlU_g Q 5 g 

2.40 

Not marginalized 


Table 2. The second column shows the current best-fit values and la uncertainties on the three- 
flavor oscillation parameters assuming normal hierarchy in the fit [52]. The third column shows 
the true values of the oscillation parameters used to simulate the ‘observed’ data set. The fourth 
column depicts the range over which sin^ 023 and 5 cp are varied while minimizing the fo obtain 
the final results. In our calculations, serves as an input parameter and then we estimate the 

value of using the relation given in Eq. (3.4) (see text for details). In the third column, we 

take i5cp = 0° to calculate the value of Amlg from Am^^. 


column of Table 1, the values of Vcc have been estimated using Eq. (3.3) where we take 
the line-averaged constant Earth matter densities® for both the baselines: p = 2.87 g/cm^ 
for the DUNE baseline and p = 3.54 g/cm^ for the LBNO baseline. Table 1 shows that 
around first oscillation maximum, the strengths of the terms , Vcc-, and 14/i are 

®These line-averaged constant Earth matter densities have been estimated using the Preliminary Refer¬ 
ence Earth Model (PREM) [115]. 
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comparable for both the set-ups under consideration even if Oe/i is as small as 10“^^. It 
means that these three terms can interfere with each other having significant impact on 
the oscillation probability which we are going to study in this section with the help of 
analytical expressions. Before we start deriving our approximate analytical expressions for 
the effective mass-squared differences and mixing angles in matter in the presence of long- 
range potential, let us take a look at the current status of the oscillation parameters. The 
second column of Table 2 shows the present best-fit values and la errors on the three-flavor 
oscillation parameters assuming normal hierarchy in the fit [52]. We use the benchmark 
values of the various oscillation parameters as given in the third column of Table 2 to draw 
the oscillation probability plots in this section and to generate the ‘observed’ data set while 
estimating the physics reach of the experimental set-ups. The ranges over which sin^ 023 and 
(5cp are marginalized while minimizing the are given in the fourth column which we will 
discuss in detail while describing the simulation method in the later section. In Table 2, 
is the effective mass-squared difference measured by the accelerator experiments 
using disappearance channel [41, 116] and it is a linear combination of Am|^ and 

The value of Am|^ is estimated from Am^^ using the relation [117, 118] 

Am|i = Am^^ -|- Am 2 i(cos^ 0 i 2 — cosdcp sin 0 i 3 sin 20 i 2 tan 023) • (3-4) 

The value of Am^^ is calculated separately for NH and IH using the above equation as¬ 
suming Am^^ = ± 2.4 X 10“^ eV^ where positive (negative) sign is for NH (IH). Note that 
through out this paper, whenever we vary 5cp or 023 or both, we calculate a new value for 
Amg]^ using Eq. (3.4). 


3.1 Analytical Expressions for the Effective Oscillation Parameters 

In a CP-conserving scenario (dcp = 0°), the effective Hamiltonian in the flavor basis given 
in Eq. (3.1) takes the form 

Hf = R23 ( 023 ) Ri3 {O13) Ri 2 {O12) H0RJ2 ( 012 ) 7 ?r 3 (^is) R^3 (^23) + k" , ( 3 . 5 ) 


where Hq = Diag (0, A 21 , A 31 ) with A 21 = Am2i/2E and A 31 = Am'^i/2E. In the above 
equation, V = Diag (Fee + 14/^,—ke/.t, 0) for Lg — symmetry. We can rewrite 77/ in 
Eq. (3.5) as 


where 


Hf = A 


31 


011 012 013 
012 022 O23 
\ Oi3 O23 O33 


Oil = A -|- IT -|- sin^ 013 + a cos^ 0 i 3 sin^ 0 i 2 , 

012 = [cos 013(a cos 012 sin 012 -|- sin 0 i 3 — osin^ 0 i 2 sin 013 )] 

v 2 


013 = [cos 013(—a cos 012 sin 012 -|- sin 0 i 3 — asin^ 0 i 2 sin 013 )] 

v 2 


(3.6) 


(3.7) 

(3.8) 

(3.9) 
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022 = 2 ~ ”*” ~ ’ 

(3.10) 

023 = 2 [cos^ ^13 ~ « cos^ 012 + OL sin^ 012 sin^ 013 ] , (3-11) 

\ 

033 = 2 [*^ 0 ®^ ^13 + acos^ 012 + a sin 013(sin 2012 + sin^ 0i2 sin 013 )] . (3.12) 


In the above equations, we introduce the terms hh, and a which are defined as 


4 = ^ VetJ. _ ‘^EVe/j, ^ 

A31 Am|]^ ’ A31 Am|]^ ’ Am|;^ ’ 


(3.13) 


and we assume that the vacuum value of 023 is 45°. Note that we have kept the terms 
of all orders in sin 0 i 3 and a which are quite important in light of the large value of 1-3 
mixing angle as was measured recently by the modern reactor experiments. Now, we need 
to diagonalize the effective Hamiltonian ^1/ in Eq. (3.6) to find the effective mass-squared 
differences and mixing angles in the presence of the Earth matter potential (Vcc) ood the 
long-range potential (14^). We can almost diagonalize Hj with the help of a unitary matrix 


U = R 23 {e^^)Rl 3 (0^3) ^12 (0^2) , 


(3.14) 


such that 

U'^HfU ~ Diag {m\^/2E, m\j2E, m\j2E) , (3.15) 

where off-diagonal terms are quite small and can be safely neglected. The lower right 2x2 
block in Eq. (3.6) gives us the angle 0^ which has the form 


tan 20 ^ = 


cos^ 013 — a cos^ 012 -k OL sin^ 0i2 sin^ 0i3 
W E OL sin 2012 sin 0 i 3 


(3.16) 


The mixing angles 9’^ and 0j^ can be obtained by subsequent diagonalizations of the (1,3) 
and ( 1 , 2 ) blocks respectively and we get the following expressions 


tan 20 ^ = 


sin 2013(1 


a sin^ 012 ) (cos 0 ^ -t- sin 0 ^) — a sin 20 i 2 cos 0 i 3 (cos 0 ^ 
\/2 (A 3 — A — W — sin^ 013 — a sin^ 0 i 2 cos^ 013 ) 


and 


sin 0 ™) 

? 

(3.17) 


tan 20 ^ = 
cos 0 ^ [sin 2013(1 


asin^ 012 ) (cos 0 ^ — sin 0 ^) -|- a sin 20 i 2 cos 0 i 3 (cos 0 ^ -|- sin 0 ^)] 
\/2 (A 2 — Ai) 


(3.18) 


where 
A 3 = 

A 2 = 


cos^ 013 -k a cos^ 012 -k a sin^ 0 i 2 sin^ 0 i 3 — IT -k 
cos^ 013 -k a cos^ 012 -k ot sin^ 0i2 sin^ 0i3 — IT — 


IE -k a sin 20 i 2 sin 0 i 3 
cos 20 ^ 

IT -k a sin 20 i 2 sin 0 i 3 
cos 20 ^ 


(3.19) 

(3.20) 
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and 


Ai — 


(A3 + A + W + sin^ 013 + a cos^ 0i3 sin^ 012) 

(A3 — A — W — sin^ 013 — a cos^ 0i3 sin^ 012) 


cos 20”* 


(3.21) 


The eigenvalues m‘f^/2E (i = 1, 2, 3) are given by the expressions 




(As + A + W + sin^ 0i3 + a cos^ 0i3 sin^ 012) 

(A3 — A — W — sin^ 013 — a cos^ 0i3 sin^ 012) 


+ - 


cos 20^ 


(3.22) 


and 


m\j2E 


A31 

~Y 


Ai + A2 


(Ai — A2) 

cos 20*” 


m\j2E 


A31 

~Y 


Ai + A2 + 


(Ai — A2) 

cos 20*” 


(3.23) 


(3.24) 





Figure 1. The variations in the effective mixing angles with the neutrino energy E in the presence 
of the Earth matter potential (Vcc) and long-range potential (14^). The left, middle, and right 
panels show the ‘running’ of 0^, 05)31 and 0^ respectively. Here, we take L = 1300 km which 
corresponds to the Fermilab-Homestake baseline and assume NH. Plots are given for three different 
choices of the effective gauge coupling 0 (the SM case), 10“®^, and 10“®^. The vacuum values 
of the oscillation parameters are taken from the third column of Table 2 and we consider dcp = 0°. 


10 




























E [GeV] 


Figure 2. The variations in the effective mass-squared differences with the neutrino energy E in 
the presence of Vcc and Left panel shows the ‘running’ of Am§ 2 ,m(— ™ 3 ,m ~ while 

right panel is for Am| 2 ,m(= ™ 2 ,m ~ Here, we take L = 1300 km and assume NH. We give 

plots for three different choices of the effective gauge coupling aefj,'. 0 (the SM case), 10“®^, and 
10“®^. The vacuum values of the oscillation parameters are taken from the third column of Table 2 
and we consider Jcp = 0°. 


With the help of Eq. (3.16), Eq. (3.17), and Eq. (3.18), we plot the ‘running’ of the 
effective mixing angles 9^, and respectively in Fig. 1 as functions of the neutrino 
energy E in the presence of Vcc and 14^. Here, we consider L = 1300 km and NH. We 
give the plots for three different choices of the effective gauge coupling 0 (the SM 
case), 10“®^, and 10“®^. The vacuum values of the oscillation parameters are taken from 
the third column of Table 2. We do the same for the effective mass-squared differences 
(see Eqs. (3.22), (3.23), and (3.24)) in Fig. 2. The extreme right panel of Fig. 1 shows 
that 9^2 approaches to 90° very rapidly with increasing E in the presence of Vcc and this 
behavior does not change much when we introduce Ve^. This is not the case for and 9^^. 
The long-range potential 14/^ affects the ‘running’ of (see extreme left panel of Fig. 1) 
significantly. As we go to higher energies, 0^ deviates from 45° and its value decreases 
very sharply depending on the strength of ae^- In case of the effect of 14/^ is quite 
opposite as compared to as can be seen from the middle panel of Fig. 1. Assuming 
NH, as we increase E, 9^ quickly reaches to 45° (resonance point) in the presence of 14/^ 
and then finally approaches toward 90° as we further increase E. These two opposite 
behaviors of 9"^ and 9^ alter the amplitudes and the locations of oscillation maxima in 
the transition probability substantially for non-zero which we discuss in the later part 
of this section. For = 10“^^ (10“®^), the resonance occurs around 4 GeV (0.6 GeV) 
for 1300 km baseline. We can obtain an analytical expression for the resonance energy 
demanding 9Y^ = 45° in Eq. (3.17). In one mass scale dominance approximation where 
Am^i can be neglected i.e. assuming a = 0, the condition for the resonance energy {Eres) 
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takes the form: 


X 3 = A + W + sin^ 013 . 


(3.25) 


Now, putting a = 0 in Eqs. (3.19) and (3.16), we get a simplified expression for A 3 which 
has the following form 


A 3 


COS^ 013 — W + + (cos^ 013)^ 



[2cos2 0i3-kk] , 


(3.26) 


since at resonance energies, the term is small compared to cos^ 0 i 3 and can be safely ne¬ 
glected. Now, comparing Eq. (3.25) and Eq. (3.26), we get a simple and compact expression 
for the resonance energy 


cos 2013 

2 V);c + 314^ 


(3.27) 


In the absence of long-range potential Ve/j,, Eq. (3.27) gives us the standard expression for 
the resonance energy in the SM framework. Eq. (3.27) suggests that for a given baseline, 
in the presence of 14^, the resonance occurs at lower energy as compared to the SM case 
and this is exactly what we observe in the middle panel of Fig. 1. The right panel of Fig. 2 
demonstrates that the effective solar mass-squared difference ^ increases with energy 

and can be comparable to the vacuum value of at higher energies in the SM case. In 

the presence of 14 /^) ^"i^i ^ increases with energy even faster and can become quite large 
at higher energies depending on the strength of Oe/^- In the SM framework, the effective 
atmospheric mass-squared difference Am|]^ ^ does not run with energy for this choice of 
baseline (see left panel of Fig. 2). But, in the presence of 14/^, the value of Am^^ ^ enhances 
a lot depending on the strength of as energy is increased. Next, to check the accuracy of 
our approximate analytical results, we compare the oscillation probabilities calculated with 
our approximate effective running mixing angles and mass-squared differences with those 
calculated numerically for the same baseline and line-averaged constant matter density 
along it. 


3.2 Demonstration of the Accuracy of the Approximation 


The neutrino oscillation probabilities in the presence of Vcc and 14 ^ can be obtained by 
replacing the vacuum expressions of the elements of the mixing matrix U and the mass- 
square differences Am?- with their effective ‘running’ values 




Uai ( 012 ^ 0 ^ 2 ,% 


5 ^23 




Am^^- 




9 9 

TTi. — 777 • 

' j,m 


(3.28) 


Incorporating the modifications due to Vcc and Ve^, the new transition probability in a 
CP-conserving scenario can be written as 

( Am? L\ 

. (3.29) 

/ 


Using Eq. (3.29), we obtain the following expressions for the appearance and disappearance 
channels [119] 




= 4U, 




ij^2 


AE 


+ AU, 




1723 sin2 


AE 
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Figure 3. ^ Ve transition probability as a function of neutrino energy E in GeV for 1300 

km (2290 km) baseline in left (right) panels. The upper panels are for the SM case without long- 
range potential. The lower panels correspond to a^fi = 10“®^. In all the panels, we compare our 
analytical expressions (dashed curves) to the exact numerical results (solid curves) for NH and IH. 
The vacuum values of the oscillation parameters are taken from the third column of Table 2 and 
we take dcp = 0°. 


- - ' ' ( 

+2 U^,sUe3Uf,2Ue2 4 sin^-4^ sin2 


AE 


AE 


+2 U^ 3 Ue 3 U^, 2 Ue 2 sin-sill ■ 


2E 


2E 


(3.30) 




= 


= ^-^U"2 


1 - C/“2 I sm 


„ Amni 

2 -^ Jj2 

AE 

.L 


1 - C/^3 I sm- 


AE 


+ 2 ( 4 = 1 .^ 

m2 m3 I 4^ 
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Figure 4. — >■ transition probability as a function of neutrino energy E in GeV for 1300 

km (2290 km) baseline in left (right) panels. The upper panels are for the SM case without long- 
range potential. The lower panels correspond to = 10“®^. In all the panels, we compare our 
analytical expressions (dashed curves) to the exact numerical results (solid curves) for NH and IH. 
The vacuum values of the oscillation parameters are taken from the third column of Table 2 and 
we take dcp = 0°. 




. ^^21,mL . Amii L 
sin-—-sin ■ 


2E 


2E 


(3.31) 


In Fig. 3, we present our approximate —)■ r'e oscillation probabilities (dashed curves) 
as a function of the neutrino energy against the exact numerical results (solid curves) 
considering L = 1300 km (left panels) and 2290 km (right panels). We give the plots 
for both NH and IH in all the panels considering line-averaged constant Earth matter 
densities for both the baselines. The upper panels are drawn for the SM case (cie/i = 0) 
where our approximate results match exactly with the numerically obtained probabilities. 
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In the lower panels, we give the probabilities considering Oe/i = 10“^^ and hnd that our 
approximate expressions work quite well in the presence of long-range potential and can 
predict almost accurate LjE patterns of the oscillation probability. In Fig. 4, we study 
the same for —)• oscillation channel and find that our approximate expressions match 
quite nicely with the numerical results. Here, we present our analytical results for NH. We 
can obtain the same for IH by changing —)• — Am|^. Following the same procedure 

and reversing the sign of V in Eq. (3.5), we can derive the analytical expressions for anti¬ 
neutrino as well. Note that in this paper, we limit our investigation to symmetry, 

though similar procedure can be adopted for Lg — Lr symmetry. 


3.3 Discussion at the Probability Level — Neutrino Case 


In this section, we discuss in detail how the long-range potential affects the full three-flavor 
neutrino oscillation probabilities in matter considering non-zero values of dcp- In Fig- 5, 
we show the exact numerical transition probability P^e as a function of the neutrino energy 
using the line-averaged constant Earth matter densities for 1300 km (left panels) and 2290 
km (right panels) baselines. We vary dcp within the range —180° to 180° and the resultant 
probability is shown as a band, with the thickness of the band reflecting the effect of dcp 
on P^e- Inside each band, the probability for dcp = 0° case is shown explicitly by the black 
dashed line. The left panels (right panels) are for 1300 km (2290 km) baseline. In each 
panel, we compare the probabilities for NH and IH with and without long-range potential. 
In the upper (lower) panels, we take ae^ = 10“®^ {oie^i = 10“^^) for the cases with long- 
range potential. We study the same for the disappearance —)• channel in Fig. 6. 

We give the similar plots for the anti-neutrino case in appendix A. 

To explain the behavior of the oscillation probabilities in Fig. 5 and Fig. 6 at the 
qualitative level, we can simplify the analytical expressions given in Eqs. (3.30) and (3.31) 
in the following fashion. The extreme right panel of Fig. 1 suggests that —)• 1 and 

cos 6y 2 0 very quickly as we increase E in the SM case or with non-zero Oe/j- So, we 
set sin 05 ^ ~ 1 and cosO"^ ~ 0 in Eqs. (3.30) and (3.31) and obtain the following simple 
expressions: 

TT? ^ Y 

^ i^e) = Sin^ 9 T 3 sin' 20^3 sin^ , (3.32) 

and 




= 


r, 00 ™T 

1 — sin^ 29 ^ sin^ 9 ^ sin^ 


'13 ‘ 


AE 

„ r, r, ArUoi 

— sin^ 29 ^ cos^ 9 ^ sin^- —= - 

AE 

Amor, 


^13 ' 


— sin'^ 9 ^ sin^ 297 ^ sin^ 


^32,m-^ 


AE 


(3.33) 


In Fig. 5, we can see that for Ue/i = 10“^^ case (upper panels), the locations of the hrst 
oscillation maxima have been shifted toward lower energies for both the baselines and also 
the amplitudes of the first oscillation maxima have been enhanced when we assume NH. 
This can be understood from the ‘running’ of 9^, 9^^ (extreme left and middle panels of 
Fig. 1), and Am‘^ 2 m, nnd right panels of Fig. 2). As we go to higher energies, 9^!^ 
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Figure 5. The transition probability P^e as a function of neutrino energy. The band reflects the 
effect of unknown Jcp- Inside each band, the probability for 5 cp = 0° case is shown by the black 
dashed line. The left panels (right panels) are for 1300 km (2290 km) baseline. In each panel, 
we compare the probabilities for NH and IH with and without long-range potential. In the upper 
(lower) panels, we take = 10“®^ (agp. = 10“®^) for the cases with long-range potential. 


increases and 9^ decreases and there is a trade-off between the terms sin^ 9^ and sin^ 20)^ 
in Eq. (3.32). Also, the value of " ^^ 21 ,m) decreases with energy as 

Am^i „j increases by substantial amount compared to „j, which shifts the location of 

the hrst oscillation maxima toward lower energies. For IH, the value of 9^ decreases fast 
with non-zero compared to the SM case, causing a depletion in the probabilities over a 
wide range of energies. In case of = 10“^^ (lower panels), there is a huge suppression in 
the probabilities at both the baselines over a wide range of energies above 1 GeV assuming 
NH. The main reason behind this large damping in the probabilities is that 9’^ approaches 
very quickly to 90° around 1 GeV or so for Oe/i = 10“®^ (see middle panel of Fig. 1) 
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Figure 6. The transition probability as a function of neutrino energy. The band reflects the 
effect of unknown Jcp- Inside each band, the probability for dcp = 0° case is shown by the black 
dashed line. The left panels (right panels) are for 1300 km (2290 km) baseline. In each panel, 
we compare the probabilities for NH and IH with and without long-range potential. In the upper 
(lower) panels, we take = 10“®^ (agp. = 10“®^) for the cases with long-range potential. 


and therefore, sin^ 20)^ —k 0, vanishing the probability amplitnde for —)• i/g oscillation 

channel. Below 1 GeV, runs toward 45° and therefore, sin^ 20^ —?■ 1, causing the 
enhancement in the probabilities. When we take IH, quickly advances to zero, causing 
a huge damping in the probabilities at all the energies. These ‘running’ behaviors of 0^, 
0 f^, and the mass-squared differences in the presence of long-range potential as discussed 
above also affect t oscillation channel (see Fig. 6) which can be explained with 
the help of Eq. (3.33). Next, we discuss how the long-range potential due to Lg — 
symmetry modifies the expected event spectra and total event rates of the DUNE and 
LBNO experiments. 
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4 Impact of Long-Range Potential at the Event Level 

We start this section with a brief description of the main experimental features of the 
DUNE and LBNO set-ups that we use in our simulation. 

4.1 Key Features of DUNE and LBNO Set-ups 

The proposed DUNE experiment [74-78, 120] in the United States with a baseline of 1300 
km from Fermilab to Homestake mine in South Dakota is planning to build a massive 35 kt 
liquid argon time projection chamber (LArTPC) as the far detector choice. This LArTPC 
will have excellent kinematic reconstruction capability for all the observed particles, re¬ 
jecting almost all of the large neutral current background. We use the detector properties 
which are given in Table 1 of Ref. [79] . As far as the neutrinos are concerned, this facil¬ 
ity will have a new, high intensity, on-axis neutrino beam, which in its initial phase, will 
operate at a proton beam power of 708 kW, with proton energy of 120 GeV, delivering 
6 X 10^® protons on target in 230 days per calendar year. In this work, we have used the 
latest fluxes being considered by the collaboration [121]. We have assumed five years of 
neutrino run and five years of anti-neutrino run to estimate the physics capabilities of this 
set-up. 

In Europe, the proposed LBNO experiment [79-82] offers an interesting possibility to 
address the fundamental unsolved issues in neutrino oscillation physics using a baseline 
of 2290 km between CERN and Pyhasalmi mine in Finland which enables us to cover a 
wide range oi L/E choices, mandatory to resolve parameter degeneracies. The Pyhasalmi 
mine will house a giant 70 kt LArTPC as a far detector which will observe the neutrinos 
produced in a conventional wide-band beam facility at CERN. The fluxes that we use in 
our simulation have been computed assuming an exposure of 1.5 x 10^^ protons on target 
in 200 days per calendar year from the SPS accelerator at 400 GeV with a beam power 
of 750 kW [122]. For LBNO also, we assume five years of neutrino run and five years of 
anti-neutrino run. We consider the same detector properties as that of DUNE. 

4.2 Event Spectrum and Rates 

In this section, we present the expected event spectra and total event rates for both the set¬ 
ups under consideration in the presence of long-range potential. We calculate the number 
of expected electron events^*^ in the f-th energy bin in the detector using the following 
expression 

Ni = —jj / dE dEA </>(E) a,^{E) R{E, Ea) P^e{E ), (4.1) 

where 0(E) is the neutrino flux, T is the total running time, Un is the number of target 
nucleons in the detector, e is the detector efficiency, and R(E, Ea) is the Gaufiian energy 
resolution function of the detector, is the neutrino interaction cross-section which has 

^°The number of positron events can be estimated using Eq. (4.1), by considering appropriate oscillation 
probability and cross-section. The same is true for events. 
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been taken from Refs. [123, 124], where the authors estimated the cross-section for water 
and isoscalar targets. In order to have LAr cross-sections, we have scaled the inclusive 
charged current cross-sections of water by a factor of 1.06 for neutrino and 0.94 for anti¬ 
neutrino [125, 126]. The quantities E and Ea are the true and reconstructed (anti-)neutrino 
energies respectively, and L is the path length. 


DUNE (35 kt) LBNO (70 kt) 



Figure 7. Expected signal and background event spectra in the appearance channel as a function 
of the reconstructed neutrino energy including the efficiency and background rejection capabilities. 
The left panel is for the DUNE (35 kt) and the right one is for the LBNO (70 kt). In each panel, 
the thick lines correspond to the SM case, whereas the thin lines are drawn assuming a^fi = 10“®^. 
In both the panels, the solid grey (blue) vertical lines display the locations of the first (second) 
oscillation maxima. We assume dcp = 0° and NH. For other oscillation parameters, the values are 
taken from the third column of Table 2. 

In our study, we consider the Vf, and Ue appearance channels, where the backgrounds 
mainly stem from the the intrinsic fe/^e contamination of the beam, the number of muon 
events which will be misidentified as electron events, and the neutral current events. In 
Fig. 7, we show the expected signal and background event spectra as a function of re¬ 
constructed neutrino energy including the efficiency and background rejection capabilities. 
The left panel shows the results for the DUNE set-up with 35 kt far detector mass. The 
right panel displays the same for the LBNO set-up with 70 kt far detector. In both the 
panels, the thick lines correspond to the SM case, whereas the thin lines are drawn assum¬ 
ing Qfe/i = 10“^^. In both the set-ups, we can clearly see a systematic downward bias in the 
reconstructed energy for the neutral current background events due to the final state neu¬ 
trino included via the migration matrices. The solid grey (blue) vertical lines display the 
locations of the first (second) oscillation maxima. The red solid histogram shows the signal 
event spectrum. Note that in the presence of long-range potential, both the signal and 
background (intrinsic contamination and misidentified muons) event spectra get mod¬ 
ified substantially. For both the baselines, we have considerable number of signal events 
around the second oscillation maximum. But, these event samples are highly contaminated 
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with the neutral current and other backgrounds at lower energies, limiting their impact. 


Channel 

DUNE (35 kt) 

LBNO (70 kt) 

Signal 

Background 

Signal 

Background 

CC 

Int+Mis-id+NC=Total 

CC 

Int+Mis-id+NC=Total 

(NH, SM) 

590 

125+29+24=178 

1228 

115+31+29=175 

P^e (NH, SM-FLRF) 

588 

123+34+24=181 

786 

112+53+29=194 

(IH, SM) 

268 

129+29+24=182 

220 

126+31+29=186 

P^e (IH, SM-bLRF) 

108 

130+33+24=187 

49 

128+50+29=207 

P^e- (NH, SM) 

116 

43+10+7=60 

117 

33+11+13=57 

P^e (NH, SM+LRF) 

44 

44+12+7=63 

22 

34+19+13=66 

P^e (IH, SM) 

210 

42+10+7=59 

484 

30+11+13=54 

Ppe (IH, SM+LRF) 

220 

41+12+7=60 

343 

29+19+13=61 


Table 3. Comparison of the total signal and background event rates in the Ve/^e appearance 
channel for DUNE (35 kt) and LBNO (70 kt) set-ups. Here ‘Int’ means intrinsic beam contamina¬ 
tion, ‘Mis-id’ means misidentified muon events, and ‘NC’ stands for neutral current. For the cases 
denoted by ‘SM-I-LRF’, we take = 10“®^. The results are shown for both NH and IH assuming 
i^CP = 0°- For both the set-ups, we assume five years of neutrino run and five years of anti-neutrino 
run. 


Channel 

DUNE (35 kt) 

LBNO (70 kt) 

Signal 

Background 

Signal 

Background 

CC 

NC 

CC 

NC 

P^^ (NH, SM) 

4889 

24 

5222 

29 

P^^ (NH, SM+LRF) 

5806 

24 

8949 

29 

P^^ (IH, SM) 

4882 

24 

5203 

29 

P^^ (IH, SM+LRF) 

5569 

24 

8519 

29 

Pp-^ (NH, SM) 

1751 

7 

1936 

13 

Ppp (NH, SM+LRF) 

2012 

7 

3257 

13 

Ppp (IH, SM) 

1752 

7 

1923 

13 

Ppp (IH, SM+LRF) 

2063 

7 

3309 

13 


Table 4. Comparison of the total signal and background event rates in the disappearance 

channel for DUNE (35 kt) and LBNO (70 kt) set-ups. For the cases denoted by ‘SM-I-LRF’, we 
take a^fj, = 10“®^. The results are shown for both NH and IH assuming ^cp = 0°- For both the 
set-ups, we consider five years of neutrino run and five years of anti-neutrino run. 


In Table 3, we show a comparison between the total signal and background event rates 
in the Ve/^e appearance channel for DUNE (35 kt) and LBNO (70 kt) set-ups. For both 
the set-ups, we assume five years of neutrino run and five years of anti-neutrino run. For 
the cases denoted by ‘SM-I-LRF’, we take = 10“^^. The results are shown for both 
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NH and IH assuming (5cp = 0°. The Earth matter effects play an important role for both 
the baselines which is evident from the fact that in the neutrino channel, the number of 
expected events is quite large compared to the IH case and in the anti-neutrino channel, 
the situation is totally opposite where we have larger event rates for IH than for NH. The 
relative difference between the number of events for NH and IH is larger for the CERN- 
Pyhasalmi baseline than the FNAL-Homestake baseline, since the impact of matter effects 
is more significant at the 2290 km baseline compared to the 1300 km baseline. In the 
presence of long-range potential with a benchmark choice of aefi = 10“®^, qualitatively, 
the trend remains the same as mentioned above. Table 3 clearly shows that in all the 
cases, the most dominant contribution to the background comes from the intrinsic t'e/^e 
beam contamination. Note that though the total signal event rate for the DUNE set-up 
in the neutrino mode with NH, does not change much due to long-range potential with 
UejM = 10“^^, but, the shape of the signal event spectrum gets affected by considerable 
amount as can be seen from Fig. 7, which enables us to place tight constraints on as 
we discuss in the results section. In our simulation, we also include the information coming 
from the disappearance channels. For these type of channels, neutral current events 

are the main source of background. Table 4 shows the total signal and background event 
rates in the disappearance channels for both the set-ups, considering hve years of 

neutrino run and five years of anti-neutrino run. For the cases marked by ‘SM-I-LRF’, 
we take = 10“®^ and we present results for both NH and IH assuming dcp = 0°- 
Interestingly, the disappearance channels are also quite sensitive to the long-range 

potential and in all the cases, we see a significant change in the total signal event rates 
with cig^ = 10“^^ as compared to the SM case. Also, the rates are different for NH and 
IH in the presence of long-range potential. The disappearance channels also play an 

important role to constrain the atmospheric oscillation parameters in the fit. 

4.3 Bi-events Plot 

In this section, we make an attempt to unravel the impact of long-range potential with 
the help of bi-events plot. In Fig. 8, we have plotted r'g vs. hg appearance events, for 
DUNE (left panel) and LBNO (right panel), considering both NH and IH and with and 
without long-range potential. Since dcp is not known, events are generated for the full 
range [—180°, 180°], leading to the ellipses. For the cases labelled by ‘SM-I-LRF’, we take 
aefjL = 10“®^. We generate these plots with sin^ 023 = 0.5 as mentioned in Table 2. The 
ellipses in Fig. 8 suggest that both the set-ups can discriminate between NH and IH at high 
conhdence level, irrespective of the choice of 5cp j and the presence of long-range potential 
with = 10“®^ does not spoil this picture. We can see from the left panel that for the 
DUNE set-up, the anti-neutrino (neutrino) event rates get reduced for NH (IH) with LRF 
as compared to the SM case. But, for CERN-Pyhasalmi baseline (see right panel) with 
more matter effect, both the neutrino and anti-neutrino event rates get diminished for NH 
and IH in the presence of long-range potential as compared to the SM case. Now, let us 
make an attempt to understand this behavior. In the presence of long-range potential, 
the locations of the first oscillation maxima shift towards lower energies (see the upper 
panels of Fig. 5 and Fig. 16), where both the fluxes and the interaction cross-sections are 
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DUNE (35 kt) 


LBNO (70 kt) 




Figure 8. Bi-events and Ue appearance) plots for NH and IH with and without long-range 
potential. The ellipses are due to all possible dcp values. For the cases labelled by ‘SM-I-LRF’, we 
take Oep = 10“®^. The left panel is for the DUNE set-up (35 kt) and the right panel is for the 
LBNO set-up (70 kt). 

small. On the contrary, at higher energies, we see a suppression in the probability where 
we have most of the neutrino fluxes and the cross-sections are also high at these energies. 
These opposite bahaviors are responsible for the large depletions in the event rates. Fig. 8 
also portrays that the asymmetries between the neutrino and anti-neutrino appearance 
events are largest for the combinations: (NH, (5cp = —90°) and (IH, <5cp = 90°). One 
striking feature emerging from both the panels is that all the ellipses get shrunk in the 
presence of long-range potential, reducing the differences in the number of events due to 
the CP-conserving and CP-violating phases. It ultimately affects the CP-coverage for the 
leptonic CP-violation searches, where we study the choices of the CP phase, dcp which can 
be distinguished from both 0° and 180° at a given confidence level. The right panel shows 
that this effect is more prominent for the LBNO set-up, severely limiting its discovery reach 
for CP-violation which we discuss in detail in the results section. 

5 Simulation Method 

In this section, we give a brief description of the numerical technique and analysis proce¬ 
dure which we adopt to estimate the physics reach of the experimental set-ups. We have 
made suitable changes in the GLoBES software [127, 128] to obtain our results. The entire 
numerical analysis is performed using the full three-flavor oscillation probabilities. Unless 
stated otherwise, we generate our simulated data considering the true values of the oscil¬ 
lation parameters given in the third column of Table 2. These choices of the oscillation 
parameters are well within their lu allowed ranges which are obtained in recent global fit 
analysis [52]. In the fit, we marginalize over test sin^ 023 and test dcp in their allowed 
ranges which are given in the fourth column of Table 2, without assuming any prior on 
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these parameters. We also marginalize over both the hierarchy choices in the fit for all 
the analyses, except for the mass hierarchy discovery studies where our goal is to exclude 
the wrong hierarchy in the fit. We keep @13 fixed in the fit as the Daya Bay experiment is 


expected to achieve a relative Icr precision of ~ 3% by the end of 2017 [129], and needless 
to say that the global oscillation data will severely constrain @13 beyond the Daya Bay 


limit before these future experiments will come online. For the atmospheric mass-squared 
splitting, we take the true value of = ± 2.4 x 10“^ eV^ where positive (negative) 

sign is for NH (IH), and we do not marginalize over this parameter in the fit since the 
projected combined data from the currently running T2K and NOz^A experiments will be 
able to improve the precision in jAm^^j to sub-percent level for maximal 623 [130]. On top 
of the standard three-flavor oscillation parameters, we have also the LRF parameter ae^ 
which enters into the oscillation probability. As far as the true and test values of Oe^ are 
concerned, we vary our choices in the range 10“®“^ to 10“®^, where the lower limit^^ corre¬ 
sponds to the cases where the oscillation probabilities almost overlap with the SM cases for 
the set-ups that we consider in this work. We take the upper limit of ae^ as 10“®^ which 
covers all the existing bounds on this parameter available from the oscillation experiments 
as discussed in section 2.2. Based on the techniques discussed in Refs. [131, 132], we use 
the following function in our statistical analysis: 


n 



(5.1) 


where n is the total number of bins and 



Above, ae^}) is the predicted number of signal events in the i-th energy bin for a 

set of oscillation parameters uj and a particular value of ae^- is the number 

of background events in the f-th bin where the charged current backgrounds are dependent 
on UJ and and the neutral current backgrounds do not depend on the oscillation 
parameters and The quantities and vr^ in Eq. (5.2) are the systematic errors on 
the signal and background respectively. For both the set-ups, we consider tt® = 5% and 
TT^ = 5% in the form of normalization error for both the appearance and disappearance 
channels. The quantities and are the “pulls” due to the systematic error on signal 
and background respectively. We incorporate the data in Eq. (5.1) through the variable 
Xi = + N^, where is the number of observed charged current signal events in the 

i-th energy bin and is the background as mentioned earlier. To estimate the total 
we add the contributions coming from all the relevant channels in a given experiment 
in the following way 



2 


+ x: 


,2 


(5.3) 


the range of the long-range force is equal or larger than onr distance from the galactic center, then 
the collective long-range potential due to all the electrons in the galaxy becomes significant. In such cases, 
these experimental set-ups can be sensitive to even lower values of Ue/j [88]. 
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Figure 9. Expected bounds on from the DUNE (35 kt) and LBNO (70 kt) experiments in 
the scenarios when the data show no signal of LRE. Results are given for four different choices of 
true values of i5cp- The left panel is for CP-conserving choices: (5cp(true) = 0° (solid lines), 180° 
(dashed lines). The right panel is for maximal CP-violating choices: i5cp(true) = 90° (solid lines), 
—90° (dashed lines). In all the cases, we assume NH as true hierarchy. 

where we assume that all these channels are completely uncorrelated, all the energy bins 
in a given channel are fully correlated, and the systematic errors on signal and background 
are fully uncorrelated. Finally, x^otai marginalized in the fit over the oscillation param¬ 
eters, both the hierarchy choices, the LRF parameter ae^ (as needed), and the systematic 
parameters as mentioned above to obtain 

6 Results 

In this section, we report our main findings. First, we present the expected constraints on 
from the proposed DUNE and LBNO experiments. Next, we quantify the discovery 
reach for of these future facilities. Then, we address how the flavor-dependent LRF, 
mediated by the extremely light Lg — gauge boson can affect the CP-violation searches 
and the mass hierarchy measurements at these upcoming facilities. 

6.1 Expected Constraints on the Effective Gauge Coupling OgAt 

In this section, we estimate the upper bounds on Og/i from the proposed DUNE and LBNO 
experiments if there is no signal of LRF in the data. This performance indicator corresponds 
to the new upper limit on Og^ if the experiment does not see a signal of LRF in oscillations. 
We simulate this situation in our analysis by generating the data at Q:g^(true) = 0 and fitting 
it with some non-zero value of Og^ by means of the x^ technique as outlined in section 5. 
The corresponding obtained after marginalizing over oscillation parameters (023, 

(5cp, and mass hierarchy) and systematic parameters in the fit, is plotted in Fig. 9 as a 
function of Ug^ (test), which gives a measure of the sensitivity reach of the DUNE or LBNO 
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Expt. 

<5cp(true) 

90% C.L. 

3ct 

Expt. 

(5cp(true) 

90% C.L. 

3cr 

DUNE 
(35 kt) 

0“ 

2.0 X 10-®®(NH) 
2.1 X 10-®®(IH) 

3.0 X 10-®®(NH) 
3.1 X 10-®®(NH) 

LBNO 
(70 kt) 

0® 

7.0 X 10-®'‘(NH) 
7.8 X 10-®'*(IH) 

1.1 X 10-®®(NH) 

1.2 X 10-®®(IH) 

180® 

2.1 X 10-®3(NH) 
2.0 X 10-®®(IH) 

3.1 X 10“®®(NH) 
3.0 X 10-5®(IH) 

180® 

7.6 X 10“®'‘(NH) 
7.0 X 10-®'‘(IH) 

1.2 X 10“®®(NH) 
1.1 X 10-®®(IH) 

90® 

1.7 X 10-®3(NH) 

1.8 X 10-®3(IH) 

2.7 X 10-®®(NH) 

2.8 X 10-®®(IH) 

90® 

7.0 X 10-®'‘(NH) 
7.0 X 10-®'*(IH) 

1.1 X 10-“(NH) 
1.1 X 10-®®(IH) 

-90® 

1.9 X 10-®3(NH) 
1.5 X 10-®®(IH) 

3.0 X 10-®®(NH) 
2.4 X 10-®®(IH) 

-90® 

7.8 X 10-®'‘(NH) 
6.0 X 10-®'*(IH) 

1.2 X 10-®®(NH) 
9.6 X 10-®'*(IH) 


Table 5. The expected bounds on from the DUNE (35 kt) and LBNO (70 kt) experiments if 
there is no signal of LRF in the data. The results are presented at 90% and Scr confidence levels 
for four different choices of true values of 5cp: 0°, 180°, 90°, and —90°. For each i5cp(true) value, 
we show the results for both the choices of true hierarchy: NH and IH. 


set-up to the effective gauge coupling of the LRF. New limits are given for four different 
choices of true values of hcp- The left panel is for CP-conserving choices: (5cp(t^^ue) = 
0° (solid lines), 180° (dashed lines). The right panel is for maximal CP-violating choices: 
(5cp(true) = 90° (solid lines), —90° (dashed lines). In all the cases, we assume NH as true 
hierarchy. Fig. 9 clearly shows that the LBNO set-up with 70 kt detector mass can place 
better limits on ae^ as compared to the DUNE set-up with 35 kt detector and the limits 
are not very sensitive to the choice of unknown (5cp(true) for both the set-ups. Table 5 lists 
the precise upper limits on which are expected from these future facilities if there is no 
trace of LRF in the data. We present the bounds at 90% (1.64 (t) and 3a confidence levels^^ 
for four different choices of true values of 5cp: 0°, 180°, 90°, and —90°. For each (5cp(true) 
value, we give the results for both NH and IH as true hierarchy choice. For an example, 
if (Icp(true) = —90° and true hierarchy is NH, then the 90% C.L. limit from the DUNE 
(LBNO) experiment is < 1.9 x 10“®^ (7.8 x 10“^^), suggesting that the constraint from 
the LBNO experiment is ~ 2.4 times better than the DUNE set-up^^. This future limit 
from the DUNE (LBNO) experiment is ~ 30 (70) times better than the existing limit^^ 
from the SK experiment [85] which is also mainly sensitive to the atmospheric mass scale 
like long-baseline experiments. At 3a, we see the same relative improvement in the LBNO 
experiment in constraining compared to the DUNE set-up. Table 5 also suggests that 
the limits are not highly dependent on the true choices of mass hierarchy. The variation 
in the upper limits on ccg^ is also not significant as we vary (Icp(true) in the range -180° 
to 180° as can be seen from Fig. 10, where we give the expected bounds at 3a and 5fT 
confidence levels from both the set-ups assuming NH as true hierarchy. Next, we discuss 


calculate this, we use the relation na — \/AyW(. In [133], it was shown that the above relation is 
valid in the frequentist method of hypothesis testing. 

^®We have checked that the larger detector mass (two times) in the LBNO set-up compared to the DUNE 
set-up is partially responsible for this improvement in the sensitivity, but also the larger path length with 
more matter effect plays an important role in this direction. 

^^This limit is quite old and was derived in a two-flavor scheme assuming = 0° [85]. One needs to 
revise this limit using the presently available full data set from the SK experiment in light of the non-zero 
and large 6 ^ 13 . 
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Figure 10. Constraints on as a function of true value of (5cp assuming NH as true hierarchy. 
Results are shown for DUNE (35 kt) and LBNO (70 kt) at Scr and 5 ct confidence levels. 


the discovery reach for ae^ if we find a positive signal of LRF in the expected event spectra 
at DUNE and LBNO. 

6.2 Discovery Reach for agfj. 

How good are our chances of observing a positive signal for LRF and hence Og/j in these 
proposed facilities? We answer this question in terms of the parameter indicator which we 
call the “discovery reach” of the experiment for We define this performance indicator 
as the expected lower limit on true values of Og^ above which the projected data at DUNE 
or LBNO would give us a signal for LRF at a certain confidence level. To find these limiting 
values, we simulate the data for various true values of Og^ and fit it with a predicted event 
spectrum corresponding to Og^ = 0. We marginalize over 623 , Sqp, mass hierarchy, and 
systematic parameters in the fit to estimate the resultant ‘^Xpiscovery which is plotted in 
Fig. 11 for DUNE (35 kt) and LBNO (70 kt) set-ups, considering four different choices of 
true values of dcp- In the left panel, we take the CP-conserving choices: (5cp(true) = 0° 
(solid lines), 180° (dashed lines). In the right panel, we consider the maximal CP-violating 
choices: (Icp(true) = 90° (solid lines), —90° (dashed lines). In all the cases, we take NH as 
true hierarchy. The nature of the curves in Fig. 11 are quite similar to the curves which 
are shown in Fig. 9, and LBNO with 70 kt detector has better discovery reach for Og^ as 
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Figure 11. showing the discovery reach for cteji expected from the DUNE (35 kt) and 

LBNO (70 kt) set-ups in the cases when one finds a signal of LRF in the data. Results are given for 
four different choices of true values of 5 cp- The left panel is for CP-conserving choices: 5cp(ti'ue) = 
0° (solid lines), 180° (dashed lines). The right panel is for maximal CP-violating choices: 5cp(ti'ue) 
= 90° (solid lines), —90° (dashed lines). In all the cases, we assume NH as true hierarchy. 


Expt. 

(5cp(ti'ue) 

90% C.L. 

3f7 

Expt. 

5cp(tme) 

90% C.L. 

3cr 


0 “ 

3.8 X 10-5“(NH) 
3.7 X 10-““(IH) 

6.5 X 10-““(NH) 
7.2 X 10-5“(IH) 


0 “ 

1.2 X 10-““(NH) 
1.4 X 10-53(IH) 

2.2 X 10-““(NH) 
2.6 X 10-““(IH) 

DUNE 
(35 kt) 

180“ 

2.9 X 10-“3(NH) 
4.0 X 10-““(IH) 

5.0 X 10-53(NH) 
7.8 X 10-5“(IH) 

LBNO 
(70 kt) 

180“ 

9.5 X 10-5^(NH) 
1.4 X 10-53(IH) 

1.9 X 10-5%NH) 
2.6 X 10-““(IH) 


90“ 

2.0 X 10-5“(NH) 
2.3 X 10-““(IH) 

5.4 X 10-®%NH) 

4.4 X 10-5“(IH) 


90“ 

9.0 X 10-®'‘{NH) 
1.0 X 10-5“(IH) 

2.0 X 10-53(NH) 
2.0 X 10-““(IH) 


-90“ 

2.0 X 10-5“(NH) 
1.8 X 10-53(IH) 

3.7 X 10-“%NH) 

4.8 X 10-®“(IH) 


-90“ 

1.0 X 10-®%NH) 
7.6 X lO-s^IH) 

2.0 X 10-53(NH) 
1.5 X 10-““(IH) 


Table 6. The discovery reach for as expected from the DUNE (35 kt) and LBNO (70 kt) 
experiments if the data show a signal of LRF. The results are presented at 90% and 3cr confidence 
levels for four different choices of true values of i5cp: 0°, 180°, 90°, and —90°. For each 5cp(true) 
value, we show the results for both the choices of true hierarchy: NH and IH. 


compared to DUNE with 35 kt, like in the case of constraints on ae^- In Table 6, we give 
the precise lower limits on true values of Oe^ which can be separated from = 0 in the 
fit at 90% and 2>a confidence levels. The results are given for both the set-ups and for four 
different choices of true values of dcp: 0°, 180°, 90°, and —90°. For each (5cp(true) value, 
we show the results for both the choices of true hierarchy: NH and IH. If we compare the 
entries in Table 6 and Table 5, then we can see that the values of the discovery reach for 
aefj, are slightly different than the constraints on Oe/i at a given confidence level and for the 
same choices of true oscillation parameters. Also, we can see that the values of discovery 
reach are marginally dependent on the choices of true dcp and mass hierarchy as we have 
seen for the constraints in the previous section. 
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Figure 12. CP-violation discovery reach as a function of true value of <5cp assuming NH as true 
hierarchy. In the left panel, we show the results for DUNE (35 kt) and the right panel is for LBNO 
(70 kt). For the ‘SM’ case, ae^ = 0 in the data and also in the fit. For each (5cp(true), we also 
give the results generating the data with three different true values of Uefi which are mentioned in 
the figure legends. In all these three cases, in the fit, we marginalize over test values of in its 
allowed range. The rest of the simulation details are exactly similar to the ‘SM’ case (see text for 
details). 


True Hierarchy 

DUNE (35 kt) 

LBNO (70 kt) 

SM 

ae/.i(t™e) = 10 

SM 

ae/iitme) = 10 

2cj C.L. 

NH (true) 

0.67 

0.62 

0.71 

0.56 

IH (true) 

0.68 

0.59 

0.73 

0.44 

3cj C.L. 

NH (true) 

0.48 

0.41 

0.55 

0.30 

IH (true) 

0.53 

0.37 

0.60 

0.12 


Table 7. Fraction of <5cp(true) for which a discovery is possible for CP-violation from DUNE (35 
kt) and LBNO (70 kt) set-ups at 2a and 3cr confidence levels. We show the coverage in 5cp(true) 
for both the choices of true hierarchy: NH and IH. For the ‘SM’ cases, we consider = 0 in the 
data and also in the fit. We also give the results generating the data with ae^(true) = 10“®^ and 
in the fit, we marginalize over test values of in its allowed range. The rest of the simulation 
details are exactly similar to the ‘SM’ case (see text for details). 


6.3 How Robust are CP-violation Searches in Presence of LRF? 

This section is devoted to study how the long-range potential due to Lf. — symmetry 
affects the CP-violation search which is the prime goal of these future facilities. Can 
we reject both the CP-conserving values of 0°, 180° at a given conhdence level? The 
performance indicator “discovery reach of leptonic CP-violation” addresses this question 
and obviously, this measurement becomes extremely difficult for the Jcp values which are 
close to 0° and 180°. In Fig. 12, we present the CP-violation discovery reach of DUNE 
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(left panel) and LBNO (right panel) as a function of true value of (5cp assuming NH as true 
hierarchy. In this plot, we generate our predicted event spectrum (data) considering the 
true value of dcp as shown in the x-axis, along with the other true values of the oscillation 
parameters given in the third column of Table 2. Then, we estimate the various theoretical 
event spectra assuming the test dcp to be the CP-conserving values 0° and 180°, and by 
varying simultaneously 023 in its 3 cj allowed range and both the choices of mass hierarchy. 
We calculate the between each set of predicted and theoretical event spectra using 
the procedure described in section 5. The smallest of all such values; Ay^py is 
plotted in Fig. 12 as a function of (lcp(true) in the range -180° to 180°. In both the panels, 
the solid red lines depict the ‘SM’ case where cte^ = 0 in the data and also in the fit. 
For each (5cp(true), we also give the results generating the data with three different true 
values of Og/x which are mentioned in the figure legends. In all these three cases, in the 
fit, we also marginalize over test values of in the range 10“^'^ to 10“^^ along with the 
other three-flavor oscillation parameters as discussed before. Fig. 12 clearly shows that 
the CP-violation discovery reach can be altered by substantial amount as compared to 
the ‘SM’ case depending on the true choice of Og^. In case of Q;g^(true) = 6 x 10“®^, we 
see a large suppression in the CP-violation discovery reach of DUNE (left panel) in the 
range 45° < (lcp(ti'ue) < 135°. We have checked that this mainly happens due to the 
marginalization over 023 in the fit where we vary sin^ 623 over a wide range (0.38 to 0.64) 
without imposing any prior on it. In case Q;g^(true) = 10“®^, the LBNO set-up (right panel) 
suffers a large depletion in the CP-violation discovery reach which can be easily explained 
with the help of bi-events plot (Fig. 8) shown in section 4.3. In the right panel of Fig. 8, 
we have seen a large reduction in the ly and i? event rates for CERN-Pyhasalmi baseline 
with Og^ = 10“®^ and this is true for both NH and IH. The differences in the number of 
events for the CP-conserving and CP-violating phases get reduced as the ellipses in Fig. 8 
get shrunk in the presence of LRF, severely deteriorating the CP-violation discovery reach 
of LBNO as can be seen from the right panel of Fig. 12. Table 7 also validates this result, 
where we compare the precise fraction of (lcp(true) for which a discovery is possible for 
CP-violation from LBNO (70 kt) and DUNE (35 kt) at 2cj and 3 ct confidence levels. For 
LBNO set-up with true NH and Q;g^(true) = 10“®^, the coverage in (icp(true) at 3cr C.L. 
reduces to 30% from 55% as we have in the ‘SM’ case. In case of true IH, the impact of 
long-range potential is even more dramatic for these future facilities. At 3a with ae^(true) 
= 10“®^, their CP-violation reach is quite minimal; only 12% for LBNO and 37% for 
DUNE while in the ‘SM’ framework, the coverage is 60% for LBNO and 53% for DUNE. 
Since, the sign of the long-range potential 14^ is opposite for neutrino and anti-neutrino, 
it affects the neutrino and anti-neutrino oscillation probabilities in different fashion. This 
feature introduces fake CP-asymmetry like the SM matter effect and severely limits the 
CP-violation search in these long-baseline facilities which can be clearly seen from Table 7. 

Finally, to see the complete picture, the fraction of (5cp(ti’ue) for which a discovery is 
possible for CP-violation is shown in Fig. 13 as a function of true value of Ofg^ assuming 
NH as true hierarchy. In each panel, we compare the performances of DUNE (35 kt) and 
LBNO (70 kt) which are shown by red and blue lines respectively. We give the results at 2fT 
(left panel) and 3a (right panel) confidence levels. In both the panels, the solid horizontal 
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Figure 13. Fraction of 5cp(true) for which a discovery is possible for CP-violation is plotted 
as a function of true value of Uep, assuming NH as true hierarchy. In each panel, we compare 
the performances of DUNE (35 kt) and LBNO (70 kt) which are shown by red and blue lines 
respectively. We give the results at 2cr (left panel) and 3 ct (right panel) confidence levels. In both 
the panels, the solid lines portray the ‘SM’ scenario where aepi = 0 in the data and also in the fit. 
For the ‘SM+LRF’ case (dashed lines), the data is generated with the true value of aefj. as shown 
in the x-axis and in the fit, we marginalize over test values of Uep, in its allowed range. The rest of 
the simulation details are exactly similar to the ‘SM’ case (see text for details). 


lines depict the ‘SM’ case where ae^ = 0 in the data and also in the fit. For the ‘SM+LRF’ 
case (dashed lines), the data is generated with the true value of Og/i as shown in the x-axis 
and in the fit, we marginalize over test values of ae^ in its allowed range. The rest of the 
simulation details are exactly similar to the ‘SM’ case as discussed before. For the values 
close to ae^(true) = 10“®^, the event spectra in the data is almost similar to the ‘SM’ case, 
but since we allow to vary in the fit in the range 10“^^ to 10“^^ along with the other 
three-flavor oscillation parameters as discussed before, we see a small suppression in the 
fraction of (5cp(true). In both the panels, around ag^(true) = 6 x 10“®^, the CP-violation 
discovery reach of DUNE deteriorates substantially, which we also observe in Fig. 12, and 
the marginalization over 023 is mainly responsible for this as we have already discussed. 
Once ag^(true) approaches toward 10“®^, the coverages in dcpltrue) for which CP-violation 
can be observed, shrink very rapidly for both the set-ups, and ultimately around ag^(true) 
= 2 X 10“^^, the coverages almost become zero. We can understand this feature from our 
discussions in section 3.1, where we have seen that in the presence of Ve/i, as we increase E, 
quickly approaches toward 45° (see middle panel of Fig. 1), and the resonance occurs at 
much lower energies as compared to the SM case. Finally, 6 ^ reaches to 90° as we further 
increase E, and the —?• Pg oscillation probabilities vanish for most of the energies where 

we have significant amount of neutrino flux. It causes a huge suppression in the event 
rates and as a result, the sensitivity goes to zero. Next, we turn our attention to the mass 
hierarchy discovery potential of DUNE and LBNO. 
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Figure 14. Discovery reach for mass hierarchy as a function of test assuming NH in the data 
and IH in the fit. The upper (lower) panels are for DUNE (LBNO). We give the results for four 
different true values of 5 cp in each panel. The solid horizontal lines in each panel show the ‘SM’ 
case where = 0 in the data and also in the fit. For the ‘SM+LRF’ case (dashed lines), the data 
is generated with the true value of = 10“®^ (left panels), 10“®^ (middle panels), and 10“^^ 
(right panels). Then, in the fit, we vary the test values of while marginalizing over 023 and i5cp- 
The rest of the simulation details are exactly similar to the ‘SM’ case (see text for details). Also, 
note that the ranges in the x-axis and y-axis are different in some of the panels. 


6.4 Impact of LRF on Mass Hierarchy Measurements 

The large Earth matter effects at both the DUNE and LBNO baselines enhance the sepa¬ 
ration between the oscillation spectra of NH and IH, and hence, we have large differences in 
the event rates for NH and IH, leading to unprecedented sensitivity toward neutrino mass 
hierarchy. Now, it would be quite interesting to see how robust are these measurements in 
the presence of LRF? A ‘discovery’ of the mass hierarchy is a discrete measurement and 
is defined as the ability to exclude any degenerate solution for the wrong (fit) hierarchy 
at a given confidence level. For hierarchy sensitivity, we first assume NH to be the true 
hierarchy and we choose a true value of dcp and Oe/x- We compute the NH event spectrum 
for these assumptions and the other true values of the oscillation parameters (see the third 
column of Table 2) and label it to be data. Then, we estimate the various theoretical event 
spectra assuming IH and a test value of as shown in the x-axis of Fig. 14, and by vary¬ 
ing simultaneously test 023 in its Scr allowed range and test <5cp in the full allowed range 
(-180° to 180°). Next, we compute the between each set of predicted and theoretical 
event spectra using the numerical technique described in section 5. The smallest of all 
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Figure 15. Discovery reach for mass hierarchy as a function of test ae^ assuming IH in the data 
and NH in the fit. The upper (lower) panels are for DUNE (LBNO). We give the results for four 
different true values of 5 cp in each panel. The solid horizontal lines in each panel show the ‘SM’ 
case where ae^ = 0 in the data and also in the fit. For the ‘SM+LRF’ case (dashed lines), the data 
is generated with the true value of = 10“®^ (left panels), 10“®^ (middle panels), and 10“^^ 
(right panels). Then, in the fit, we vary the test values of while marginalizing over 023 and i5cp- 
The rest of the simulation details are exactly similar to the ‘SM’ case (see text for details). Also, 
note that the ranges in the x-axis and y-axis are different in some of the panels. 


such values: ^Xmb shown in Fig. 14 as a function of Q;e/. 4 (test) for given choices 
of ae/i(ti'ue) and (5cp(ti'ue). As mentioned above, we always assume NH in the data and 
IH in the fit while generating the curves in Fig. 14. The upper panels portray the per¬ 
formance of DUNE (35 kt), while the lower panels are for LBNO (70 kt). In each panel, 
the resnlts are given for four different choices of (5cp(true) and the solid horizontal lines 
depict the ‘SM’ case where is zero in the data and also in the fit. For the ‘SM-I-LRF’ 
case (dashed lines), the data is generated with the true value of as mentioned in the 
top part of each panel, and in the fit, we vary the test values of Og^ while marginalizing 
over 023 and Jcp- The rest of the simulation details are exactly similar to the ‘SM’ case 
as mentioned above. Though in Fig. 14, we have shown the results for three benchmark 
values of ag^(true) = 10“®^ (left panels), 10“^^ (middle panels), and 10“®^ (right panels), 
but, we have checked that for DUNE, the mass hierarchy sensitivity always stays above the 
standard expectations irrespective of (5cp(t3'3ie) provided the true value of Og^ < 5 x 10“^^. 
There is a large suppression in the appearance event rates when we generate the data with 
a true value of Og^ around 5 x 10“^^, and if we further increase the valne of ag^(true). 
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the statistical strength of the data reduces very rapidly, and the sensitivity goes below the 
standard expectation. The upper right panel in Fig. 14 clearly shows this behavior. In case 
of LBNO, the mass hierarchy discovery reach never goes below the ‘SM’ value irrespective 
of (5cp(true) if the true choice of is smaller than 10“^^. Once we consider the true value 
of cteij, > 10 “®^, the appearance event rates get reduced in data by considerable amount, 
causing a significant drop in the sensitivity which can be clearly seen from the lower middle 
and right panels in Fig. 14. In Fig. 15, we generate the data with IH and fit it with NH. 
We see almost similar behavior in all the panels of Fig. 15 as we have noticed in Fig. 14. 


7 Summary and Conclusions 

Flavor-dependent long-range leptonic forces mediated by the extremely light and neutral 
bosons associated with gauged Lg — or Lg — symmetries, constitute a minimal exten¬ 
sion of the SM preserving its renormalizability and can lead to interesting phenomenological 
consequences. For an example, the electrons inside the Sun can generate a flavor-dependent 
long-range potential I4^/gr at the Earth surface which can give rise to non-trivial three 
neutrino mixing affects in terrestrial experiments, and could influence the neutrino propa¬ 
gation through matter. The sign of this potential is opposite for anti-neutrinos, and affects 
the neutrino and anti-neutrino oscillation probabilities in different fashion. This feature 
invokes fake CP-asymmetry like the SM matter effect and can severely affect the leptonic 
CP-violation searches in long-baseline experiments. In this paper for the first time, we 
have investigated in detail the possible impacts of these long-range flavor-diagonal neutral 
current interactions in the oscillations of neutrinos and anti-neutrinos in the context of 
future high-precision superbeam facilities, DUNE and LBNO. The key point here is that 
for long-baseline neutrinos, Am^/2E ~ 2.5 x 10“^^ eV (assuming Am^ ~ 2.5 x 10“^ eV^ 
and E ~ 5 GeV) which is comparable to 14/^ even for ofg^ ~ 10“®^, and can influence 
the long-baseline experiments significantly. For the Fermilab-Homestake (1300 km) and 
CERN-Pyhasalmi (2290 km) baselines, the Earth matter potentials are also around 10“^^ 
eV (see Table 1 ), suggesting that Vcc can also interfere with Ug^ and Am|]^/ 2 E, having 
substantial impact on the oscillation probability. We have explored these interesting pos¬ 
sibilities in detail in this work. In this paper, we have presented all the results considering 
the Lg — L^ symmetry. Similar analysis can be performed for the Lg — L^ symmetry which 
we will present elsewhere. 

We have derived approximate analytical expressions for the effective neutrino oscilla¬ 
tion parameters to study how they ‘run’ as functions of the neutrino energy in the presence 
of both long-range and Earth matter potentials. We have also obtained a compact and 
simple expression for the resonance energy, where 0 i 3 becomes 45° in the presence of both 
Vcc and Vefj,. We have observed that in the presence of Ug^, as we increase the neutrino 
energy, ^13 in matter quickly approaches toward 45°, and the resonance occurs at much 
lower energies as compared to the SM case. Finally, ^13 in matter reaches to 90° as we 
further increase the energy, causing a large suppression in the appearance probability for 
most of the energies where we have significant amount of neutrino flux for both the set-ups. 
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As a result, the event rates get reduced which can be clearly seen from the bi-events plot 
in Fig. 8. 

As the long-range potential due to gauged Lg — symmetry can change the standard 
oscillation picture of these future facilities significantly, we can expect to place strong 
constraints on if these experiments do not observe a signal of LRF in oscillations. For 
an example, if 5cp(true) is —90° and true hierarchy is NH, then the expected bound from 
the DUNE (35 kt) set-up at 90% C.L. is Ogfi < 1.9 x 10“^^. The same from the LBNO (70 
kt) experiment is Ugfi < 7.8 x 10“®^, suggesting that the constraint from LBNO is 2.4 times 
better than DUNE. This future limit from the DUNE (LBNO) experiment is almost 30 
(70) times better than the existing bound from the SK experiment [85]. We have noticed 
that these future limits on ag^ from DUNE and LBNO are not very sensitive to the true 
choice of (5cp and mass hierarchy. We have also estimated the discovery reach for agf^ if we 
find a positive signal of LRE in the expected event spectra at DUNE and LBNO. We have 
found that the spectral information on the signal and background events is quite crucial 
to constrain/discover this new long-range force. 

We have also studied in detail the CP-violation discovery reach of DUNE (35 kt) and 
LBNO (70 kt) in the presence of LRE. We have seen that the CP-violation measurements 
can be deteriorated by considerable amount as compared to the standard expectation 
depending on the true value of Ug^. At 3a with ae;x(true) = 10“®^ and true NH, the 
coverage in (5cp(true) for which a discovery is possible for CP-violation is 41% (30%) for 
DUNE (LBNO) while in the standard case, the coverage is 48% for DUNE and 55% for 
LBNO. In case of true IH, the impact of long-range potential is even more striking for 
these future facilities. As an example, if aeAi(true) = 10“®^, their chances of establishing 
CP-violation are quite minimal: only for 37% (12%) values of (5cp(ti'ue), DUNE (LBNO) 
can reject both the CP-conserving values 0° and 180° in the fit at 3a, while in the ‘SM’ 
framework, DUNE (LBNO) can do so for 53% (60%) values of true dcp- As the true value 
of Oe/i approaches toward 10“®^, the coverages in (lcp(true) for which CP-violation can be 
observed, diminish very quickly for both the set-ups, and ultimately around ae^(true) = 
2 X 10“®^, the coverages almost become zero. 

Pinally, we have asked the question, how robust are mass hierarchy measurements in 
these future facilities in the presence of LRE? In the standard case, due to the large Earth 
matter effects at the Eermilab-Homestake and CERN-Pyhasalmi baselines, both DUNE 
and LBNO can resolve the issue of mass hierarchy at very high confidence level. Now, if 
LRE exists in Nature, then for DUNE, the mass hierarchy sensitivity remains above the 
standard expectations provided the true value of < 5 x 10“^^. In case of LBNO, the 
mass hierarchy discovery reach does not go below the ‘SM’ value as long as the true value 
of agf_i is smaller than 10“®^. 
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A Discussion at the Probability Level — Anti-Neutrino Case 
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Figure 16. The transition probability P^g as a function of anti-neutrino energy. The band reflects 
the effect of unknown dcp- Inside each band, the probability for (5cp = 0° case is shown by the 
black dashed line. The left panels (right panels) are for 1300 km (2290 km) baseline. In each panel, 
we compare the probabilities for NH and IH with and without long-range potential. In the upper 
(lower) panels, we take = 10“®^ (agp. = 10“®^) for the cases with long-range potential. 

In Fig. 16, we plot the exact numerical transition probability 12^ —)• z2g as a function of 
anti-neutrino energy. The band shows the impact of unknown (5cp- Inside each band, 
the probability for <5cp = 0° case is shown by the black dashed line. The left panels 
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(right panels) are drawn for 1300 km (2290 km) baseline. In each panel, we compare the 
probabilities for NH and IH with and without long-range potential. In the upper (lower) 
panels, we consider = 10“®^ (ag^ = 10“®^) for the cases with long-range potential. 
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Figure 17. The transition probability Ppp as a function of anti-neutrino energy. The band reflects 
the effect of unknown i5cp. Inside each band, the probability for (5cp = 0° case is shown by the 
black dashed line. The left panels (right panels) are for 1300 km (2290 km) baseline. In each panel, 
we compare the probabilities for NH and IH with and without long-range potential. In the upper 
(lower) panels, we take = 10“®^ (og/i = 10“®^) for the cases with long-range potential. 


Fig. 17 shows the exact numerical —)• disappearance probability as a function of 
anti-neutrino energy. The thin band portrays the mild impact of unknown <5cp • Inside each 
band, the probability for (5cp = 0° case is given by the black dashed line. The left panels 
(right panels) are drawn for 1300 km (2290 km) baseline. In each panel, we compare the 
probabilities for NH and IH with and without long-range potential. In the upper (lower) 
panels, we consider = 10“®^ = 10“®^) for the cases with long-range potential. 
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